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The phenomenon of periodic orbit scarring of eigenstates of 
classically chaotic systems is attracting increasing attention. 
Scarring is one of the most important "corrections" to the 
ideal random eigenstates suggested by random matrix the- 
ory. This paper discusses measures of scars and in so doing 
also tries to clarify the concepts and effects of eigenfunction 
scarring. We propose a new, universal scar measure which 
takes into account an entire periodic orbit and the linearized 
dynamics in its vicinity. This measure is tuned to pick out 
those structures which are induced in quantum eigenstates 
by unstable periodic orbits and their manifolds. It gives en- 
hanced scarring strength as measured by eigenstate overlaps 
and inverse participation ratios, especially for longer orbits. 
We also discuss off-resonance scars which appear naturally on 
either side of an unstable periodic orbit. 



I. INTRODUCTION 
A. Background 

The modern field of quantum chaology often associates 
classically chaotic motion on the one hand with aspects of 
random matrix theory (RMT) on the other jlj . These as- 
pects include level repulsion in the quantum spectrum as 
given by the appropriate random matrix ensemble, Gaus- 
sian random wavefunctions with Bessel correlations, etc. 
Since Hamiltonian dynamics cannot be truly random, nu- 
merous recent contributions to the field address the many 
sorts of "corrections" to the random matrix approxima- 
tion. One of those corrections is the phenomenon of scar- 
ring of quantum eigenstates by isolated unstable periodic 
orbits of the corresponding classical system |E| . 

In the early 1980's MacDonald in unpublished work ||] 
found clear evidence of non-isolated marginally unstable 
periodic orbits in certain stadium eigenstates (which he 
named the "bouncing ball" states). He also tentatively 
noted the possible influence of an isolated, unstable pe- 
riodic orbit on a few of the calculated eigenstates, but 
gave no further attention to this effect. In the subsequent 
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first published account of numerically computed stadium 
eigenstates by MacDonald and Kaufman Q (which to 
the authors' knowledge contains the first eigenstates re- 
ported for any completely chaotic system), attention was 
focussed not on periodic orbit effects but on the nodal 
structure of the eigenstates. The conclusion was that the 
expected pattern of nodal lines for random eigenstates 
had been reached. This early work also noted broad 
agreement with the Bessel function form of the eigenstate 
spatial self-correlation function for chaotic billiards, as 
first discussed and predicted by Berry || . The paper ]4j 
appeared during the first epoch of quantum chaos the- 
ory, when much excitement was being generated by not- 
ing similarities between computed eigenstate properties 
of classically chaotic Hamiltonian systems and RMT [0 . 

Stronger numerical evidence for the influence of in- 
dividual periodic orbits on eigenstates, together with a 
theoretical explanation for scarring in a chaotic system, 
was introduced in 1984 j2j. Loosely speaking, a scar is a 
concentration of extra and unexpected (as compared to 
the RMT prediction) eigenstate density near an unsta- 
ble classical periodic orbit. This extra concentration has 
no classical analogue, which puts scarring into the family 
of quantum localization effects. A semiclassical theory 
for the existence and strength of scars was given, using 
time domain arguments and dynamics linearized around 
the periodic orbits. This theory has seen a number of 
extensions and applications, including Bogomolny's co- 
ordinate space theory of scarring || and Berry's Wigner 
phase space theory 0. (The 1984 paper [|| had been es- 
sentially a Husimi phase space, or Gaussian wavepacket, 
theory.) All these theories were based on the linearized 
dynamics in the vicinity of a periodic orbit, but there 
were important differences. For example, an essential 
ingredient to the strength of scarring, the Lyapunov sta- 
bility exponent of the periodic orbit, enters only in the 
wavepacket approach |^], while the important observa- 
tion of "knots" of high density at self-conjugate (focal) 
points in coordinate space along the orbit was made by 
Bogomolny 

Scarring has been shown to affect physical systems of 
various sorts ^ and even the performance of devices such 
as a tunnel diode H . Recently one of us showed that the 
decay of mctastable states can be strongly affected by 
scarring, in that highly anomalous lifetime distributions 
are possible depending on where decay channels are lo- 
cated with respect to the shortest periodic orbits of the 
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system JTo[ . 

The following litany of properties of eigenfunction scar- 
ring has led to some confusion and to several attempts 
at developing quantitative measures of this phenomenon: 
(1) The manifestations of scarring can be subtle or ob- 
vious. (2) Measures of scarring can be basis dependent. 
(3) Scar-like structures are found to occur even in artifi- 
cially constructed purely random wavefunctions [fll| . (4) 
Statistical fluctuations allowed by RMT might account 
for some apparently scarred states. 

Recent work in our laboratory fll2|-[i4j has focussed 
much attention on quantifying the phenomenon of scar- 
ring, confirming the role of the instability exponent of 
a given periodic orbit, and further examining the conse- 
quences of these findings, including experimental issues 
and the effects of antiscars. Related work of ours has 
noted that scarring is concordant with bounds on the er- 
godicity of eigenstates as developed by Schnirelman [fL5| , 
Zelditch |ll| , and Colin de Verdiere |L7) . A review of re- 
cent developments in the theory of scarring was recently 
written by one of the authors Q]. The present work 
completes an important part of the picture by explicitly 
addressing the basis dependence of measures of scarring 
and arriving at a universal and optimal basis for mea- 
suring scars, while not diminishing the utility of simpler 
and more ad hoc measures. 



B. Measures 

If eigenstates were ideal random matrix states, then 
all probe states would be equivalent, in that Gaussian 
random statistics in one basis guarantees Gaussian ran- 
dom statistics in every other. In the RMT literature, 
it is sometimes noted that for any given member of the 
ensemble there will be a diagonalizing basis; however, 
this basis is non-generic and is itself randomly varying 
from one member of the ensemble to another. One way 
of approaching the corrections to RMT in Hamiltonian 
systems with a classical analogue is to show that there ex- 
ist special bases which are nonrandom and which come 
from deterministic dynamical evolution. These special 
bases bring the Hamiltonian into a manifestly non-RMT 
form. Any basis which systematically shows non-RMT 
wavefunction statistics for a classically chaotic system is 
thus potentially significant. Seen in this light, the ba- 
sis dependence of scar measures should be expected and 
even exploited. 

One such special basis is that of complex Gaussian 
wavepackets. Complex Gaussians have adjustable posi- 
tion and momentum expectation values, and satisfy min- 
imum uncertainty conditions in some system of axes in 
phase space, making them excellent measuring devices 
for the structure of eigenstates in phase space. This basis 
was the one chosen in [|| . For Gaussians centered on pe- 
riodic orbits, asymptotically exact (h — > 0) semiclassical 
dynamics for a fixed short time places rigorous non-RMT 



constraints on the statistics of eigenstate projections onto 
the Gaussian. Some of the eigenstates (precisely which 
ones cannot be specified in this short-time theory) are 
then required to have large projections onto the periodic 
orbit centered Gaussian (these are the scarred states), 
while many more are shown to have anomalously small 
projections (the anti-scarred states). We review why this 
is so in the next section. The inverse participation ratio 
(IPR) of such orbit-centered Gaussian packets is anoma- 
lous, and is governed by the classical Lyapunov stability 
exponent A as 

IPR~i (1) 

for small A. (Note that fi does not appear in this scaling, 
implying survival of the scarring phenomenon into the 
classical limit.) Some individual eigenstate projections 
onto the Gaussian basis were shown to be enhanced by 
at least j- over the RMT expectation (|] (again, exactly 
which ones is not known and would require much longer 
time information). It is sometimes stated that scar the- 
ory is not a theory of individual eigenstates. While that 
is true in many respects, especially of the energy averaged 
approaches such as Bogomolny's the Husimi (Gaus- 
sian packet) phase space theory of scars j|) predicts there 
must exist individual scarred states, especially for small 
instability exponent A. 

As versatile as the Gaussian basis is, there are choices 
to be made and certain optimizations possible which fur- 
ther sensitize the probe basis to the structures which clas- 
sical dynamics imprints onto the eigenstates of classically 
chaotic systems. Before discussing this further, we note 
a reductio ad absurdum which places restrictions on how 
far the refinement of scarring measures can go. Some 
years ago, Tomsovic and Heller [jL8j were successful in 
constructing a high quality scarred eigenstate of the sta- 
dium billiard using only semiclassical methods (the over- 
lap of the semiclassical state with the exact eigenstate 
being 0.87). Using such near-eigenstates as a probe ba- 
sis would lead to extreme non-RMT behavior in which 
all but one eigenstate have small projection onto the test 
state. Furthermore, including large parts of the classical 
invariant manifolds leading far from the region of a given 
periodic orbit subverts the idea of a scar of a periodic 
orbit. 

Fortunately there is a quite natural stopping point in 
the construction of a test basis: we use only the linearized 
dynamics (tangent map) near any given periodic orbit in 
constructing measures of scarring. In this way we arrive 
at test states that are understandable in terms of sim- 
ple invariant manifold structures of classical phase space 
near periodic orbits. Although the test states can be 
more complicated than a single Gaussian, these more so- 
phisticated test states are still determined by short time 
linear dynamics (of order of the time over which the dy- 
namics is linearizable, which scales as | \ogh\/X). (Recall 
that a given orbit can be highly nonlinear, yet possess a 
linearizable tangent map in its vicinity.) 
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In constructing linear scar measures we still have a 
number of choices to make: 

• What is the uncertainty zone of the test state in 
phase space (in the case of a Gaussian, this is the 
uncertainty ellipse in phase space). Scar measures 
will change for example depending on whether the 
ellipse lies along the stable or unstable manifolds of 
the orbit or away from them. 

• For a given period-one periodic orbit, do we con- 
struct a test state with a Gaussian placed at one 
point along the orbit or do we construct a tube 
(closely spaced Gaussians on the orbit) of some 
sort, and is this tube to be a coherent superpo- 
sition of the Gaussians or an incoherent one? (In 
the case of a discrete-time map, this corresponds 
to placing a Gaussian at each periodic point of a 
periodic orbit with length greater than one.) 

• Should we take coherent linear combinations of a 
given Gaussian and its pre- and post-images, pro- 
ducing a new test state consisting of several differ- 
ent Gaussians at each point along the orbit? (This 
idea leads to the "universal" measure of scarring.) 



C. Brief history of scar measures 

There are several threads in the attempt to make good 
measures of scarring. The original approach Q amounted 
to projection onto single Gaussians (the Husimi mea- 
sure); an 0(I/A) enhancement in the infinite-time aver- 
age return probability for a Gaussian placed on an unsta- 
ble periodic orbit was noted for small A. This implies that 
some eigenstates in specified energy ranges are systemat- 
ically enhanced by 0(I/A) in the periodic orbit regions 
over the RMT predictions. Much later it was realized 
that this local enhancement has a dramatic effect on the 
tails of the |f/'| 2 distribution 13|. 



Any theory of scarring implies some measure of the 
effect. The first theory of wavefunction scarring in po- 
sition space was developed by Bogomolny ||. Bogo- 
molny smoothed the wavefunction intensity over some 
small energy range AE using the semiclassical Green's 
function; scars are represented as smoothed sums over 
effectively finitely many periodic trajectories of the sys- 
tem. Bogomolny's semiclassical Green's function ap- 
proach is very closely related to our wavepacket dynamics 
method, as the semiclassical Green's function can be ob- 
tained from the semiclassical time-domain propagator by 
a stationary-phase Fourier transform. One difference be- 
tween the approaches is that Bogomolny envisions sum- 
ming over a large number of periodic orbits to get as close 
as possible to an energy domain resolution of order of a 
mean level spacing. As mentioned above, in some sys- 
tems it is indeed possible to use semiclassical methods to 
compute individual eigenstates of the system (T^] . In fact 
for this purpose one needs information only about orbits 



of period up to the mixing time (which scales logarith- 
mically with K) rather than the Heisenberg time (which 
scales as a power law). However, our aim here is to make 
predictions about the distribution of scarring strengths 
based only on linearized information around one periodic 
orbit; for this purpose most other orbits which produce 
additional oscillations in the density of states may be 
treated statistically |20|. It is important to note in this 
context that if we are measuring wavefunction intensities 
on a given short classical periodic orbit T 3 , then in the 
semiclassical limit there are no other short orbits that 
come close to this orbit (on a scale of K) in phase space. 
The only oscillatory contributions which will need to be 
taken into account are from orbits closely related to or- 
bits homoclinic to V [homoclinic orbits are those that 
approach V at large negative times, perform an excur- 
sion away from V into other regions of phase space, and 
then again approach V at large positive times]. In fact, 
in the fi — * limit the periodic orbit sum for a point x 
on a given periodic orbit V can be written equivalcntly 
as a contribution from the orbit V itself plus a sum over 
trajectories homoclinic to V. Although the two points of 
view (periodic and homoclinic sum) are mathematically 
equivalent, the homoclinic sum approach makes explicit 
the special role of the orbit V near which we are mak- 
ing measurements. In the homoclinic return formalism, 
it is also straightforward to see that the long-time re- 
currences of a wavepacket launched at x are correlated 
and enhanced in a way that is determined entirely by the 
stability matrix of the short orbit V '. 

A position space basis, though obviously physically 
natural in many measurement situations, is not gener- 
ally an optimal one for detecting scar effects. Unless the 
periodic point x happens also to be a focusing point of 
classical trajectories near the orbit, only a small fraction 
of the total scar strength is captured in the position ba- 
sis, and the fraction becomes smaller as h decreases (or 
as the energy increases). An easy way to see this is to no- 
tice that the effects of a classical trajectory in quantum 
mechanics generically extend to a region around the or- 
bit scaling not as a wavelength but rather as the square 
root of a wavelength (and similarly the affected region 
scales as the square root of the total number of chan- 
nels in momentum space). Thus, unless either the stable 
or unstable manifold of the orbit V at periodic point x 
happens to be oriented along the momentum direction, 
the position space basis will not be optimal, as reflected 
in the falling off of the focusing prefactor with energy in 
the semiclassical Green's function (and similarly the mo- 
mentum basis will not be optimal, unless one of the two 
invariant manifolds is oriented along the position direc- 
tion). All this will become more clear in the exposition 
of the following section. In any case, one should keep in 
mind that a position space basis can always be considered 
as a special limiting case of the Gaussian wavepacket test 
state, where the position uncertainty of the wavepacket 
becomes comparable to a wavelength, and the momen- 
tum uncertainty becomes large. 
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A Wigner phase space analysis of the scarring phe- 
nomenon was given by Berry |7[. Berry considered the 
Wigner function, again smoothed over an energy inter- 
val AE near E. Being formulated in phase space, the 
approach more closely resembles that of . Working in 
Wigner phase space instead of Husimi space also elimi- 
nates the need to choose the (apparently arbitrary) ec- 
centricity and orientation of the Gaussian wavepackets. 
The downside of Wigner phase space is the absence of 
a positivity condition on the Wigner distribution; thus 
the value of the spectral function cannot be considered 
as corresponding to an intensity or a probability of being 
found near a certain point x (and random matrix the- 
ory is therefore not applicable). The Husimi function, 
which is manifestly positive definite, is identically a phase 
space smoothing of the Wigner distribution over a phase 
space region scaling as ft. The ambiguity in choosing the 
Gaussian centered on x over which this smoothing is to 
be performed is indeed an important issue, to be con- 
sidered carefully in the following. We will see that to 
obtain the maximal scarring effect, the Gaussian must 
be chosen to be properly oriented along the stable and 
unstable directions at the periodic point. [An arbitrarily 
large wavepacket width is allowed along either of these 
directions, with a correspondingly small width in the or- 
thogonal direction. Also, strong, but non-maximal, scar- 
ring will generally be obtained for any wavepacket with 
width scaling as y/h in both the position and momentum 
directions.] 

A common limitation of the analyses |§,|]j7) is that they 
make no prediction about the properties of the spectral 
fluctuations on scales much smaller than H/Tp, where 
Tn ~ Tp I A is the decay time of the unstable orbit with 
period Tp. Therefore it is not possible to make quanti- 
tative predictions about specific individual wavefunction 
intensities, participation ratios, etc., without explicitly 
doing a Gutzwiller sum over all periodic orbits. Even if 
the sum can be performed, it is by no means clear that it 
will converge in all cases (e.g. in systems where caustics 
are important [f3). When the sum does converge it may 
produce individual semiclassical wavefunctions very dif- 
ferent from the quantum eigenstates, due to diffraction 
and other "hard quantum" effects. Furthermore, such 
Heisenberg-time calculations are extremely sensitive to 
small perturbations on the system. What one would like 
is to be able to say precisely how often a given single- 
wavefunction scar strength will appear on a given orbit, 
at what energy, and at what parameter values. In the 
semiclassical limit, this can in fact be done using only in- 
formation about linearized dynamics near the orbit itself, 
and, in some cases, about a few strong isolated homo- 
clinic recurrences which cannot be treated statistically. 

Agam and Fishman j^] define the weight of a scar 
by integrating the Wigner function over a narrow tube 
in phase space, of cross-section fi, surrounding the peri- 
odic orbit. Li and Hu integrate over coordinate space 
tubes |23|. Alternatively, de Polavieja, Borondo, and 
Benito p4[ construct a test state highly localized on a 



given periodic orbit using short-time quantum dynam- 
ics. 

Klakow and Smilansky |2^| have used a scattering ap- 
proach to quantization to study the wavefunctions of bil- 
liard systems. They treat carefully the wavefunctions on 
the Poincare surface of section, and relate their properties 
to scarring in configuration space. Ozorio de Almeida ]26| ] 
uses the Weyl representation to establish connections 
between classical and quantum dynamics, with particu- 
lar application to the semiclassical Wigner function and 
scars. Tomsovic [ p7| has used parametric variation as a 
new method for studying scar effects; scars are shown to 
induce correlations between wavefunction intensities on a 
periodic orbit and the level velocities of these wavefunc- 
tions when certain system parameters are varied. We also 
mention the work of Arranz, Borondo, and Benito p§| ] 
who have probed the intermediate region between regular 
and strongly chaotic quantum behavior, and have shown 
how scarred states first arise from the mixing of pairs of 
regular wavefunctions as h is decreased (but well before 
one reaches the semiclassical limit which is the main focus 
of the present work) . Finally, several groups |^,[3lJ have 
studied the hyperbolic scar structures associated not only 
with the periodic orbit itself but with its invariant man- 
ifolds and homoclinic orbits. 

In the next section we discuss scarring as measured by 
individual Gaussian wavepackets, which was the basis of 
A single localized test state may be optimized to 
conform to the classical invariant manifolds in the vicin- 
ity of a scar. In subsequent sections we go considerably 
beyond this measure, refin ing our templates to better 
detect scarring. In Section III we address the apparent 
arbitrariness in the choice of a point along the orbit at 
which to make the measurement, and in the eccentric- 
ity of the test Gaussian, and eliminate these ambiguities 
by building a wavepacket-averaged measure of scarring. 
Following this, in Section |v| , we use coherent linear com- 
binations of the localized test states as a more sensitive 
measure. In Section 0, off-resonance scars living on ei- 
ther side of an unstable periodic orbit are shown to follow 
naturally from our formalism. In Section VI extensions to 
higher-period orbits and continuous time are discussed, 
followed by concluding remarks in Section VII 



II. GAUSSIAN WAVEPACKET SCARRING 

A. Semiclassical dynamics of a Gaussian wavepacket 

We begin with a review of the original (Gaussian 
wavepacket) theory of scarring, as discussed in detail re- 
cently in |L2]]. In the course of the discussion, the key 
concepts of the autocorrelation function, the short-time 
spectral envelope, nonlinear recurrences, and the inverse 
participation ratio will be introduced. We will also see 
the inherent limitations of measuring scar strength using 
single Gaussian test states, pointing the way to the con- 
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struction of improved "scarmometers" [] in the following 
sections. 

Consider an arbitrary (unstable) periodic orbit of a 
chaotic system. For the purpose of simplifying the expo- 
sition, and without loss of generality, we take the periodic 
orbit to be a fixed point of a discrete-time area-preserving 
map on a two-dimensional phase space. If the periodic 
orbit in question is in fact a higher-period orbit of such 
a map, or is an orbit of a continuous-time dynamics in 
two spatial dimensions, we can reduce the problem to the 
preceding case by iterating the original map, or by taking 
a surface of section map, respectively^] 

We start with a fixed point at the origin of phase space. 
Furthermore, we can take the stable and unstable direc- 
tions at the fixed point to be vertical (p) and horizontal 
(q), respectively (we can always get the local dynamics 
into this form by first performing a canonical transfor- 
mation on the coordinates). Then the only parameter 
describing the local (linearized) dynamics near the or- 
bit is A, the instability exponent for one iteration of the 
orbit. Locally, the equations of motion are given by 



q 
p 



q' = e xt q 



P 



(2) 



We now turn to the construction of a test state which 
can be used to measure the intensity of eigenstates near 
the chosen periodic orbit. An obvious choice is a Gaus- 
sian wavepacket centered on the fixed point: 



1 



na 2 h 



1/4 



j-q 2 /2a 2 h 



(3) 



This is a minimum-uncertainty state centered at the ori- 
gin of phase space, with width a\fh in the q— direction 
and y/E/cr in the p— direction, a is at this stage an arbi- 
trary parameter: a 1 is the aspect ratio of the phase-space 
Gaussian, typically chosen to be of order unity. Ambigu- 
ity in the choice of a is an important issue that we will 
return to at the beginning of the next section. 

Eigenstate overlaps with our test state will provide 
a good measure of eigenstate intensities near the peri- 
odic orbit; however, we find it useful to begin by working 
in the time domain (our results will then be applied to 
eigenstate properties in the following subsection). 

For small enough ft, the wavepacket \a a ) and its short- 
time iterates are contained well within the linear regime. 
As long as the wavepacket stays in the phase space region 
surrounding the periodic orbit in which the linearized 
equations of motion Eq. || apply, the evolution of the 
wavepacket is completely semiclassical, given simply by 



The authors thank Eugene Bogomolny for coining this 
term. 

2 The issue of higher-period orbits and continuous time will 
be addressed explicitly in Sections [II and [vi|. 



the stretching of the q— width parameter a. More explic- 
itly, at short times we have 



U*\a a 



lin 



(4) 



where U is the unitary operator implementing the quan- 
tum discrete-time dynamics, Uu n represents the quanti- 
zation of the linearized behavior near the periodic orbit, 
and t is time, measured in units of a single mapping. 
Here —(f) is a phase associated with one iteration of the 
periodic orbit: it is given by the classical action in units 
of h, plus Maslov indices if appropriate. 

The autocorrelation function of the wavepacket is de- 
fined as the overlap of the evolved wavepacket with itself: 



A(t) = (a^lc/'K), 

which at short times is seen from Eqs. ||, || to be 

„—i<t>t 



^/cosh(Ai) 1 



(5) 



(6) 



by performing a simple Gaussian integration. The 'lin' 
subscript indicates that Eq. || describes the piece of the 
autocorrelation function coming from the linearized dy- 
namics around the periodic orbit. For a weakly unstable 
orbit (small A), An n (t) is slowly decaying, with strong 
recurrences happening for the first 0(1/A) iterations of 
the orbit. We note that the short-time autocorrelation 
function Au n (t) is cr-independent, a fact that will prove 
important later on. 

At longer times, namely beyond the log time, which 
scales as 



Ti og 



log fN 
A 



(7) 



the wavepacket leaves the linearizable region and nonlin- 
ear recurrences begin to dominate the return probability. 
Here N is the total number of Planck-sized cells in the 
accessible phase space (also equal to the dimension of the 
effective quantum mechanical Hilbert space), and / is the 
fraction of this phase space (typically O(l)) in which the 
linearized equations of motion (Eq. |^) apply. The nonlin- 
ear recurrences correspond to a piece of the wavepacket 
leaving the linear regime along the unstable manifold, 
undergoing complicated dynamics far from the periodic 
orbit, and eventually coming back along the stable man- 
ifold to intersect the original wavepacket. Semiclassi- 
cally, these recurrences are given by a sum over points 
homoclinic to the original periodic orbit (i.e. points that 
approach the periodic orbit both as i ^ +oo and as 
t — > — oo). 

Because the long-time homoclinic orbits come back 
with complicated accumulated phases, and the number 
of these recurrences grows exponentially with time, one 
might expect the total long-time return amplitudes to be 
given by Gaussian random variables. In fact, however, 
contributions from all homoclinic points lying on a single 
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homoclinic orbit (i.e. those that are exactly time-iterates 
of one another) come back in phase with each other, giv- 
ing rise to short-time correlations in A(t) for large t [ ^2[ . 
These correlations are related to the short-time dynamics 
of the original Gaussian wavepacket. In fact, we can write 
the return amplitude at times T\ og <t< Tjj {Th = N is 
the Heisenberg time, where individual eigenstates begin 
to be resolved) as a convolution 



A(t) =Y, A md{r)A lin {t- 



(8) 



Here -Au n is the short-time return amplitude, and A rn d 
has the statistical properties of an uncorrelated random 
Gaussian variable. In effect, random recurrences due to 
Gaussian fluctuations must have "echoes" that mirror the 
initial short time decay, since the recurrences re-load the 
initial state. Also, 



obtained by combining the general expression Eq. |Tl| 
with the short-time overlap dynamics of the Gaussian 
wavepacket (Eq. ^ , is fundamentally an interference phe- 
nomenon, and signals a kind of quantum localization, as 
we shall see next. Note that in the limit A — * (weak 
instability) we have 



(13) 



i.e. the enhancement factor in the long-time return prob- 
ability is proportional to A -1 



B. Local density of states 

We now define S(E) to be the fourier transform of the 
autocorrelation function, 



< A ind (T) > = 
< A* ind (T)A ind ( T ') > = —5 TT , 



(9) 



The prefactor \/N provides the proper classical normal- 
ization: in the absence of interference effects, the prob- 
ability to come back is equal to the probability for vis- 
iting any other state in the Hilbert space. The average 
in Eq. |^ is taken over long times r, T\ og <C r <C Th, 
and/or over an ensemble of systems which all have the 
same linearized dynamics around our chosen periodic or- 
bit. In either case, the total size of the Hilbert space N 
(= l/h for a phase space area normalized to unity) has 
been assumed to be large. We then obtain 

< A(t) > = 

< A*(t)A(t + A) > = 1 ^n(s)Aun(s + A) • (10) 

s 

At times beyond the Heisenberg time, this gets modified 
H to 



< A*(t)A(t + A) >= -^AU^AUs - 



-A). (11) 



F is a factor associated with the discreteness of the eigen- 
states: it is 3 for real eigenstate-test state overlaps and 
2 for complex overlaps. 

The long-time autocorrelation function is thus self- 
correlated on a scale A ~ A -1 . Qualitatively, this can be 
understood on a purely classical level: once probability 
happens to come back to the vicinity of a weakly unsta- 
ble periodic orbit, it tends to stay around before leaving 
again. On the other hand, the overall enhancement in 
the total return probability at long times: 



< l^)| 2 >= 



F 



OO 

N ^ 

s= — oo 



1 



cosh(As) ' 



(12) 



S(E) 



1 

27 



E 



i. El 



(14) 



For a non-degenerate spectrum, it is easy to see (by in- 
serting complete sets of eigenstates) that 



S( J E)=^|{n|a,)| 2 ^- J B„), 



(15) 



where E n are the eigenvalues of the dynamics, and |n) 
are the corresponding eigenstates. Thus, we obtain the 
local density of states at the wavepacket \a a ) by fourier 
transforming its autocorrelation function A(t). Cutting 
off the sum in Eq. 14| at ±Ti og , or equivalently by includ- 
ing only linearized dynamics around the periodic orbit, 
we obtain the smoothed local density of states: 



Sun(E) 



E 

t 



A Vin (t)e 



iEl 



(16) 



an envelope centered at quasienergy E — <p (see Eq. |(]), 
of width 8E ~ A, and of height ~ A -1 (a factor of 2-7T 
has been inserted into the definition of Sn n for future 
convenience). E = <fi is the analogue of the EBK quanti- 
zation condition for integrable systems; here, because of 
the instability of the orbit, scarred states can live in an 
energy range of O(A) around the optimal energy. States 
with energy more than 0(A| log A|) away from resonance 
tend to be antiscarred (i.e. they have less than expected 
intensity at the periodic orbit). 

Now long-time (nonlinear) recurrences as in Eq. || lead 
to fluctuations under the short-time envelope in the full 
spectrum S(E). Because these recurrences involve a ran- 
dom variable convoluted with the short time dynamics, in 
the energy domain we obtain random fluctuations mul- 
tiplying the short-time envelope. (It is easy to see phys- 
ically that the random oscillations must multiply the 
smooth envelope: if they were merely added to it, the 
total spectrum would become negative away from the 



G 



peak of the envelope.) Finally, at the Heisenberg time 
Tjj = N, individual states are resolved Jl2|,[l3|], and we 
see a line spectrum with a height distribution given by 



Ina a = \ {n\a a )\ 2 = r an S\i n (E n ) , 



(17) 



where r an are random variables (with mean < r an >= 
1 /N) drawn from a chi-squared distribution of one degree 
of freedom (two degrees of freedom for complex (n|a ff )). 
Thus, in the end we obtain a random (Porter-Thomas) 
line spectrum S(E), all multiplying the original linear 
envelope. 

Before concluding this review, we mention the notion 
of an inverse participation ratio (IPR), a very useful mea- 
sure for studying deviations from quantum ergodicity. 
We define 



(18) 



(Note that 52 ^na a = 1 by normalization.) Being the 
first non-trivial moment of the eigenstate intensity (J n a„ ) 
distribution, the IPR provides a convenient one-number 
measure of the strength of scarring (or any other kind of 
deviation from quantum ergodicity). The IPR would be 
unity for a wavepacket that had equal overlaps with all 
the eigenstates of the system; the maximum value of N 
is reached in the opposite (completely localized regime), 
when the wavepacket is itself a single eigenstate. Random 
matrix theory predicts an IPR of F, the strong quantum 
ergodicity factor defined above in Eq. [Pi] . 
From Eqs. Ill, we see that 



AT T ~ 1 
t=0 



(19) 



as one might expect, localization is associated with an 
enhanced return probability at long times. Now from 
Eq. [l^ we see that scar theory predicts an enhancement 
in the IPR over random matrix theory: 



IPR a<r =F^ 



1 



cosh(As) 



F 



A " 



(20) 
(21) 



where in the last line the limit of small A has been taken. 
(F, as before, is 3 or 2, depending on whether the states 
are real or complex, respectively.) The IPR thus decom- 
poses into a product of two contributions: the shape of 
the short-time envelope coming from the linear dynam- 
ics around the periodic orbit, and a quantum fluctuation 
factor F, as predicted by Porter-Thomas statistics. 



C. Limitations 

The analysis of the previous two subsections has been 
extensively tested in numerical studies liaO], which 



show that the statistical properties of eigenstate overlaps 
with Gaussian wavepackets can indeed be described by 
the scar theory. However, there are inherent limitations 
in this approach. An obvious one is the ambiguity in 
the choice of wavepacket width a. A wavepacket of any 
width can be used (as long as it and its short-time iter- 
ates are well-contained in the linearizable region, which 
condition is always satisfied for small enough h), result- 
ing in the same short-time overlaps, and thus in identical 
smoothed spectra 5u n . The IPR is also expected to be 
enhanced by the same factor for each such wavepacket, 
depending only on the decay exponent of the periodic 
orbit itself. It seems intuitively clear that a better mea- 
sure of scarring should be obtainable by appropriately 
combining information from wavepackets of all different 
aspect ratios, thus looking at a hyperbolic phase-space 
region surrounding the stable and unstable manifolds of 
the periodic point. Such a test state for measuring scars 
would incorporate knowledge of the full linearized dy- 
namics around the periodic point, not just knowledge 
about the location of the periodic point itself. 

The ambiguity and apparent arbitrariness of the pre- 
ceding definition of scarring seems even more pronounced 
in the case of a higher-period orbit of a map, or for a pe- 
riodic orbit of a continuous-time dynamics. In either of 
those cases, the analysis above can be performed at any 
periodic point lying on the orbit. Yet it is known from ex- 
perience that scars tend to live not at one periodic point 
only but along the entire orbit. Here, also, more infor- 
mation could presumably be gained by looking at the 
behavior of an eigenstate near all points on a periodic 
orbit instead of one only, thus obtaining a fuller measure 
of wavefunction scarring. 

The preceding ambiguities will, in the following sec- 
tions, point us towards a universal measure of wavefunc- 
tion scarring (in the regime of small A, where the effect is 
expected to be significant), a measure which takes full ad- 
vantage of the entire periodic orbit and the full linearized 
dynamics in the vicinity of this orbit. In the process, we 
will see how the insights of various earlier contributors 
to this field MM 



24| , P7|j29| can be incorporated into the 
resulting general approach. 



III. WAVEPACKET INTENSITY AVERAGING 

A. Density matrix test states 

Consider again the fixed point of a classical area- 
preserving map, as introduced in Eq. 0. Given the appar- 
ent arbitrariness in the choice of wavepacket which came 
out of our discussion in the previous section, it seems 
natural to extend our measure of scarring, replacing the 
pure Gaussian test state with a density matrix which 
gives weight to Gaussians of all widths: 



p = M \ die t2/T « \ a<ye xt){a ae xt\ 



(22) 
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We choose the exponent A to be the instability exponent 
of the unstable orbit (Eq. |^), but this choice is in fact ar- 
bitrary, and A can be reabsorbed into the overall normal- 
ization factor J\f and the time cutoff Tq. In the absence 
of the cutoff To, the hyperbolic test state would be com- 
pletely scale invariant, giving equal weight to Gaussians 
of all aspect ratios, from tall and thin to short and wide. 
The cutoff is, however, necessary because the linearized 
dynamics of Eq. [2] is in fact only valid in a finite classical 
region around the periodic orbit (and it also eliminates 
possible normalization difficulties). We note that Eq. |2^ 
is an incoherent superposition of wavepackets of differ- 
ent width, designed to give an unbiased measure of a 
scar (where the arbitrariness of a given choice of er is 
removed) . 

Let the classical region in which the dynamics is lin- 
earizable be given by a square in phase space, with area 
A (the exact shape and area are not important, as we will 
soon see). Then choosing a = 1, we see that the evolved 
width aVhe xt reaches the edge of the linearizable region 
at time 



1 , A 

To = 2A log 4?r 



(23) 



(Up to constants, To is the same as the log time discussed 
above in Eq. (r|) We see that a factor of order one am- 
biguity in the definition of the region A will lead only to 
an additive ambiguity in the cutoff time To, irrelevant in 
the semiclassical limit A 3> Ti. Of course the condition 
that the Gaussian just touching the boundary should be 
suppressed by a factor of I/e (as opposed to 1/e 2 , 1/e 3 , 
etc.) is still somewhat arbitrary, leading us to the more 
general form 



(24) 



where c is now an arbitrary constant of order unity. In 
any case, the ambiguity we previously had in the choice 
of a (which could be anywhere from \ffi to 1/y/h, a huge 
range in the semiclassical regime fiCl) has now been re- 
duced to a factor of order one constant c in the definition 
of p. 

In the following section, where we consider coherent 
superpositions of Gaussian wavepackets, we may wish to 
use a more stringent criterion for the cutoff time To, tak- 
ing into account the form of the leading nonlinearity of 
the dynamics near the periodic orbit. Thus, consider the 
worst-case scenario, where the stable and unstable man- 
ifolds both curve quadratically as we move away from 
the periodic orbit. Then in order for the curvature to 
be unimportant as the Gaussian stretches along one of 
the two manifolds, we may require that the distance by 
which the unstable manifold deviates from the horizontal 
line at position q be less than the vertical width (mo- 
mentum uncertainty) of a state with position width q. 

This means 0(q 2 ) < Sp ~ h/q, so the maximum dis- 

1/3 

tance q for which this holds scales as h ' . Assuming 



the same situation obtains along the stable manifold, 
we obtain that the linearizable area scales as h 2/3 , i.e. 
A = AtfW* for some classically selected area Aq. Then 
\og(A/h) = I \og(A /h), i.e. we may take the size of the 
linearizable region to be the ?i— independent value Aq, 
provided we also replace c — > c/3. Because of the 0(1) 
ambiguities already present in the choice of c, we will not 
dwell here on the numerical values appropriate to various 
systems. In any case, as we will discuss when subject- 
ing our results to numerical tests in the following section, 
there is always a tradeoff between larger c leading to more 
localization and smaller c giving more precise agreement 
with analytical predictions. 

Coming back now to Eq. [2^, we notice that we could 
instead have chosen a hard cutoff for the sum over Gaus- 
sians, e.g. 



/To 
r 
-Jo 



dt \a ae \t)(a 



(25) 



This would not qualitatively affect our discussion cither 
here or in the following sections (particularly Section IV , 
where we discuss coherent superpositions of Gaussian test 
states) . The form of Eq. ^2|, and its extensions which will 
follow in future sections, is however convenient because 
it allows for relatively straightforward analytical calcula- 
tions. 

Another important property of Eq. |2^ is that p has the 
same form in momentum space as in configuration space, 
as can be seen easily by noting that the momentum width 



u p — u q — u 
variant under t 



and that Eq. |22| is manifestly in- 
t. Thus, the stable and unstable 
manifolds of the hyperbolic point are treated symmetri- 
cally in our definition. 

The idea of averaging overlap intensities can of course 
be extended to resolve another difficulty we encountered 
at the end of the previous section, namely the apparent 
ambiguity in treating periodic orbits that are not fixed 
points. For a period Tp orbit of a map, we write 



P = 



Tp-l 

■E 

p=0 



\ a x p ,a) ( a x p ,<j\ j 



(26) 



where \a Xp ,a) is a wavepacket of width a along the unsta- 
ble manifold, but centered at periodic point x p instead of 
at the origin. Similarly, for a continuous-time dynamics, 
we can write 



dx 



(27) 



where the x coordinate parametrizes the periodic orbit in 
phase space, and at each periodic point the wavepacket is 
chosen to have width a x along the direction of the orbit 
and width a in the unstable direction at that point on 
the orbit: 

a x ,a,aA x ' iV') ~ ex P[ - ( X> - x ) 2 - y' 2 / " 2 ^ 

+ ip x (x' -x) /%]. (28) 
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Here {x,p x ) are the position of a phase space point on the 
periodic orbit and the corresponding momentum, while 
y' is a coordinate along the unstable manifold of the orbit 
at point (x,p x ). Eq. |2J] has been written down already in 
Rcf . , and a connection was made there to the phase- 
space tubes of Agam, Fishman et al. P2j. 

Now the orbit averaging of Eqs. 26, 27] can of course be 
combined with the width averaging introduced in Eq. |22] : 
thus, in the case of a map we may write 



We can construct a wavepacket-averaged local density of 
states analogous to Eq. |l5| 



S P (E)=Y J In P 5{E-E n ). 



and a corresponding inverse participation ratio 



(32) 



(33) 



Tp-1 



P = N^2 / dte- t2/T « \a Xpae x t }(a Xp ^ e xt\ 
p=o J 



(29) 



Here one may legitimately ask why we perform sep- 
arately the averaging at each periodic point: we could 
instead have made use of the orbit dynamics and ob- 
tained Gaussian wavepackets centered at each of the pe- 
riodic points starting with one wavepacket only and al- 
lowing it to evolve according to the linearized laws of 
motion. More explicitly, if we take \a) to be a Gaussian 
wavepacket of width a along the unstable direction, cen- 
tered at periodic point xq , we may construct a dynamical 
density matrix 



Pdyn 



+oo 



-t 2 /TlT ( 



° |ai in (*))(aii„(t)| 



(30) 



Here |ai; n (f)) is the original Gaussian evolved in accor- 
dance with the linearized dynamics: for example, if t is 
an integer multiple of the period Tp, then \au n (t)} is cen- 
tered at the same periodic point as |q), but with width 
Xt / Tp . T is defined as before (Eq. using the full 



ere 



instability exponent A for one iteration of the entire prim- 
itive orbit. X/P is the exponent per time step; hence the 



30 



factor of Tp in Eq. 

For small A, where not much stretching has taken place 
over one period of the orbit, not much difference exists 
between the averaging methods of Eq. ^6] and Eq. |3(]. We 
return to this connection between (linearized) dynamics 
and improved test states in Section IV, where coherent 



superpositions of Gaussian wavepackets are discussed. 
Here it suffices to note that because all wavepackets be- 
ing averaged over in pa yn are exact time-evolutes of one 
another (at least in the linear approximation), they all 
have exactly the same local density of states and inverse 
participation ratio. In fact this absence of real averaging 
is there even for the full matrix p of Eqs. |26|, in 
the limit of small A, as we shall see next. 



B. Measures of scarring from incoherent averaging 

Our measure for the strength of scarring for a given 
eigenstate \n) is now simply 



I np = (n\p\n) 



(31) 



Notice that S P (E) is nothing other than a weighted sum 
of the densities S(E) of Eq. and thus follows the 
same linear envelope S\i n which we have discussed in the 
previous section. The only thing possibly different about 
Sp(E) are the oscillations under this envelope. To under- 
stand how these oscillations in the averaged local density 
of states IPRp might differ from the Porter-Thomas fluc- 
tuations one finds for a single wavepacket, we need to 
study correlations between local densities of states for 
different wavepackets centered on the same periodic or- 
bit. 

In general, given two wavepackets \a) and \b), we can 
define a long-time averaged transport probability P a b fH] 

as 



T-l 



For a nondegenerate spectrum we easily see 



Pab = 



\{a\n)\ 2 \(b\n)\ 2 = Y,InaIn b - 



(34) 



(35) 



In particular, the IPR as defined in the preceding section 
corresponds to the special case \a) = \b): 

lPR a = NP aa =Nj2lL- (36) 



The P a b can be thought of as the covariance matrix of 
the densities of states for different wavepackets, with P aa 
being the variances or diagonal matrix elements; the cor- 
relation between two densities of states is then given by 



Cab — 



Pab 



VKaPl 



bb 



(37) 



We begin with the simplest case, where wavepackets \a) 
and \b) are exact time iterates of one another: \b) = \a(t)) 
for some time t. There, of course P a b = P aa — Pbb, and 
the correlation is unity (the two local densities of states 
S(E) are identical). More explicitly (from Eq. 20) we 
have in this case 



Pab = Paa = Pbb = ^ff(A) , 



(38) 



where F/N is the RMT prediction for the quantum long- 
time return probability and 
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s(A) = £|<a|LC|a>| 2 = £ 



1 



cosh(As) 



(39) 



is a scarring IPR enhancement factor (see Eq. g0|). In 
the last equality the periodic orbit in question has been 
taken to be period one (a fixed point). 

Now consider the opposite extreme case, where 
wavepackets |o) and \b) lie on different periodic orbits 
of the same classical action (and not related by any sym- 
metry) . Then the two local densities of states S a (E) and 
Sb{E) share the same linear envelope (coming from short 
time dynamics), but have completely uncorrelated long- 
time fluctuations: 



Ina — ^an5lin(-^n) 
Inb = r bnS\i n (E n ) 



(40) 



with r an , r on uncorrelated chi-squared variables with 
mean 1/N (see Eq. Fjh. Then 



Pn, 



Pbb = §5(A) 
1 



Cab 



1 

F 



(41) 



The correlation in this case is of order unity but still less 
than one. 

Finally, choose the two wavepackets \a) and \b) lying 
on the same orbit but not exact time-iterates of one an- 
other. We can think of wavepacket |6) as having a part 
composed of \a) and its short-time iterates and another 
part which is statistically independent of \a) although it 
lies on the same orbit. The fraction of |6) which is corre- 
lated with | a) is given by a (normalized) sum of squared 
overlaps of \b) with \a) and its iterates: 



Es\(b\ULW)\ 2 _ 9ab(\) 



Iin 



3(A) 



(42) 



where Eq. |42| serves as the definition of g a b(X). In both 
numerator and denominator linearized evolution is used, 



9abW = 



cosh(A(s + z)) 



(43) 



Cab = 



1 - 



F 



g Q fc(A) 
.9(A) 



1 

F 



(45) 



Now as the exponent A becomes small, any optimally 
oriented wavepacket \b) lying on the periodic orbit be- 
gins to look more and more like an iterate of any other 
wavepacket \a). In that limit, 



hm — — — 

a^o g(X) 



1. 



(46) 



and so the correlation C a b — ► 1. So the key result is 
that when scarring is strong (A — > 0), the spectra of all 
optimal wavepackets centered on the periodic orbit in 
question are virtually identical, making unnecessary any 
averaging over width or position along the orbit: 



(n\p\n) -> \(a a \n)\ 2 , 



(47) 



for an arbitrary \a a ) along the periodic orbit. Therefore, 
in this limit any wavepacket individually provides a uni- 
versal measure of scarring intensities, obviating the need 
to construct tubes and other averaging devices. 

We now proceed to examine quantitatively the behav- 
ior of C a b hi Eq. ^5|. First, however, we will introduce 
a model ensemble of systems which will allow us to test 
numerically this quantitative prediction and others ob- 
tained in the following sections. 



C. Ensemble averaging over hard chaotic systems 

The classical area-preserving map we will use for our 
'numerical experiments' is defined on the unit square 
(q,p) € [0,1] x [0,1], and consists of two parts. The 



first step is a three-strip generalized baker's map pi 13 



with strip widths wq + w\ + W2 = 1 • Each vertical strip 
i of width Wi < 1 and height 1 is stretched horizontally 
by a factor oil/ Wi and compressed vertically by a factor 
of Wi to make it into a horizontal strip of height Wi and 
width 1. The three strips are then stacked on top of each 
other (left becoming bottom and right becoming top) to 

Wi to 



reconstruct the unit square. Defining s, 
be the left edge of strip i, we have 

q' = (q~ Si)/wi 



= E 



j<i 



P = Si 



pWi 



(48) 



where |6) is related to some exact time- iterate of \a) by 
a stretch of e Xz of the Gaussian along the unstable man- 
ifold. 

We then have a sum of two contributions: 



Pab = Jff(A) 



9afc(A) 



5(A) 



^(A) 



1 



ffab(A) 

5(A) 



F - 1 1 

— ^— 5ab(A) + -^gW, 



(44) 



giving 



where the initial position q lies in the i-th strip, i.e. s, < 
q < Si+i. The second and final step is a kicked map |35| ] 
implemented in the left and right strips of the square, 
leaving the middle strip undisturbed: 



p"= p >-v.'(L 



! mod 1 



q" = Si> + [(</ - Si>) +p"wi> modi 



(49) 



Here i' denotes the number of the strip (0 or 2) containing 
q'. The entire mapping Eqs. Ea, Efl is now iterated. 
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The convenience of this two-step model lies in the fact 
that for any choice of kick potentials Vq and V2 act- 
ing on the left and right strips, respectively, the mid- 
dle strip experiences only baker-like horizontal stretch- 
ing and vertical shrinking. Thus, there is always a fixed 
point of the system in the middle strip, with coordinates 
If =Pf= W o+w2 ' and stretching exponent A = | logwi|. 
Furthermore, the stable and unstable manifolds of this 
fixed point are always locally vertical and horizontal, re- 
spectively, consistent with our canonical form Eq. ^. The 
kicked maps acting on the left and right strips serve to 
provide parameters which can be easily varied to pro- 
duce ensemble averaging over the details of the nonlinear 
long-time recurrences without affecting the local dynam- 
ics around the periodic orbit which is being studied]^] We 
choose kick potentials 



are statistical, and do not reflect finite-size effects. The 



V ,2(X) = 



1 2 

-r 



0,2 



(2*r) s 



sin 2ttx . 



(50) 



with Kq 2 arbitrary parameters. The condition |-Ko,2| < 1 
is sufficient to ensure hard chaos, without regular regions 
fl9f . More general kicking potentials could of course have 
been used, but we find that the two parameters -Ko,2 
provide a sufficiently large ensemble for our purposes. 

The system is quantized in a straightforward and con- 
ventional way, by multiplying the unitary matrices imple- 
menting baker's map and kicked map dynamics pl] , |33[ ]. 



D. Numerical tests 

We proceed to test the density of states correlations 
C a b (Eq. [45]) for the fixed point orbit of the map intro- 
duced above. The wavepacket \a) with horizontal width 
ay/K and vertical width \ffija is placed on the fixed 
point. We then define a family of wavepackets \b(z)) 
of (horizontal) widths o~e Xz y/h. Notice that for integer 
z, \b{z)) is an exact iterate of \a) (in the linear approx- 
imation), and thus in that approximation the densities 
of states are identical and the correlation C a b = 1. The 
differences are expected to be greatest at half-integer z 
where \b) is most unlike any iterate of \a). The correla- 
tion C a b is now plotted as a function of z for < z < 1 
in Fig. |l|. Two sets of data are given, differing in 
the stability exponent of the periodic orbit (which is 
easily adjusted by varying the middle strip width wi). 
The two values used were A = log 5 (upper curve) and 
A = log 10 (lower curve). In each case, the numerical 
data comes from an ensemble average over systems of 
size N = 1/h = 200. The errorbars shown in the figure 



3 For a billiard system, the analogous procedure would be to 
take a given short periodic orbit and then create an ensemble 
of systems by deforming the boundary in such a way that the 
original orbit is unaffected. 



theoretical curves are obtained from Eq. 45 and require 
only the single parameter A. The agreement between the- 
ory and data is quite good; furthermore we see just how 
large the correlations are even for not very small expo- 
nents. Thus, for an orbit with a stretching factor of 10 
per iteration (A = log 10), the correlation does not go be- 
low 0.95 even for the maximally unrelated wavepackets 

(*=y2). 

In Fig. g we plot this minimum correlation C a b(z = 
1/2) on the vertical axis, versus the scarring enhance- 
ment factor g(X) on the horizontal axis. Four data points 
are used, corresponding (from left to right) to stretching 
factors e A = 20, 10, 5, 2.5. Both the expected aver- 
age enhancement g and the inter-wavepacket correlation 
Cab{z = 1/2) are uniquely given theoretically as func- 
tions of the exponent A. Again, the data agrees very 
well with the theoretical predictions. We see that a scar- 
ring enhancement factor of 2 (corresponding to stretching 
exponent A w log 5) is associated with a minimum cor- 
relation of 0.99 between the least correlated wavepackets 
on that orbit. Strong scarring thus automatically elimi- 
nates the ambiguity in measuring the strength of Gaus- 
sian wavepacket scarring. 

The question then becomes whether it is possible in 
any way to take advantage of our knowledge of the or- 
bit and its invariant manifolds to produce a scarring 
test state that would do better than a single Gaussian 
wavepacket. Indeed, this is possible, and what is nec- 
essary is to use coherent quantum superpositions of test 
states [[50 24 27J] instead of the density matrix approach 
investigated in this section. 



IV. COHERENT WAVEPACKET SUMS: 
ENHANCED SCARRING 

A. Theory 

As suggested already in we can construct a "lin- 
earized eigenstate" as a normalized coherent sum of 
Gaussian wavepackets centered on a periodic orbit. For 
a fixed point orbit, we write 



* = N / dte- t2/T ° \a ae x t ) 



(51) 



(see Eq. |2^). To is a linearized dynamics time cutoff 
as defined in Eq. and the normalization constant M 
ensures (\I/|\I/) = 1. Just as was done for the density 
matrix in the preceding section, Eq. ^l] can be general- 
ized in a straightforward way to higher-period orbits and 
to continuous time. However, to make the presentation 
more transparent the examples here and in the follow- 
ing section are restricted to the case of a fixed point, the 



generalizations being left to Section VI 



If the dynamics away from the periodic orbit were ex- 
actly linear, we could take the cutoff Tq to infinity and ob- 
tain a stationary state with quasi-energy <j) (phase e" 1 *, 



11 



see Eq. [4j): hence the name "linearized eigenstate". In 
reality, a finite cutoff is necessary because the ratio of 
the size of the linearizable region of phase space A to % 
is finite. However, if this ratio is large (as it will always 
be in the semiclassical limit % — * 0), most of the state '3/ 
maps to itself under Eq. |J, producing a large autocorre- 
lation function at short times. In the case of Gaussian 
wavepacket scarring, the extent to which the short-time 
return probability differs from unity (and thus the extent 
to which perfect localization fails to be achieved) is deter- 
mined by the instability of the orbit (i.e. by the amount 
by which A is different form zero). As measured using 
the improved test state the absence of complete local- 
ization is given by the failure of the linearized dynamics 
at long times. 

The test state "J/ lives not only at the periodic point, 
but also along the invariant manifolds. Its autocorrela- 
tion (^^(t)) decays only on the order of the log-time 
Tq ~ A -1 log A/%, as we show explicitly below in Eq. |5^. 
This makes ^ a much sharper measure of the scar char- 
acter of an eigenstate, and for small % we expect to see 
much stronger localization as measured by \& than by an 
individual wavepacket \a a ). We will now proceed to show 
this explicitly. The construction of W is extremely sim- 
ple, requiring only one piece of information beyond what 
we already needed for the single wavepacket, namely the 
(approximate) size of the region in which the dynamics is 
linearizable. No knowledge of long-time dynamics, non- 
linear recurrences, or any other periodic orbits is needed. 

We also notice that in the strong scarring limit A — > 0, 
we could just as well have used only the linearized time 
iterates of \a) (rather than wavepackets of all widths) as 
in Eq. bl] to construct the hyperbolic test state: 



dyn 



t 



E 



,-t 2 /^J4>t 



a\in{t)} 



(52) 



(cf. Eq. |3^) . This form makes manifest the close connec- 
tion between the construction of the scarring test state 
and the linearized classical dynamics (Eq. |^). It also 
makes almost trivial the generalization to higher period 
orbits and to continuous time (see Eq. ||(] and also the 
fuller discussion in Section VI). The main disadvantage 
of the form pdyn (as compared to p) is that the former re- 
quires the arbitrary choice of initial width a. However, as 
we have seen in the previous section, this choice of start- 
ing wavepacket has no effect on any measured quantities 
in the A — > limit (where replacement of the integral by 
a sum is justified). 

We begin as in Section |n| by evaluating the short-time 
autocorrelation function 

<(t) = <*|LMt)|*>. (53) 
A straightforward calculation using Eqs. |^, [3] gives: 



■ /</>/ 



dy ■ 



v / coshy 



(54) 



The overall normalization constant Q can be fixed by 
requiring A* n (0) = = 1. The integration variable 

y is a time variable scaled by A to make it dimensionless. 
In the limit To A 3> 1, i.e. log 4 3> 1, the exponential 
simplifies and we obtain 



AUt) = Qe- m J 



dy 



e -t 2 /T* 

\/coshy 



,-i<t>t -t 2 /T* 



(55) 



Now we see explicitly that the decay rate of our test state 
'J is indeed given by the log-time scale To. Of course for 
the linear autocorrelation function to be a good measure 
of the total return amplitude A (i), even for times t less 
than To, the state *f> must be well contained inside the 
linear region. This can be done by adjusting the constant 
c in Eq. ^4|: numerically we will see below that good 
quantitative agreement with Eq. |54| is obtained for c ~ 
0.6. In any case, the precise value of this constant does 
not affect any of the important scaling arguments which 
will follow. 

As in Section ||, the inverse participation ratio has 
an enhancement factor associated with the short time 



recurrences: 



IPR* 



<I(H*)I 4 > 

< |(n|*)| 2 > 2 

Fj2\AUt)\ 2 
t 



(56) 
(57) 
(58) 



where in the last line the limiting form Eq. ^5] has 
been used, and To taken to be large. In the Gaussian 
wavepacket case, the IPR scaled with the orbit decay 
time A -1 ; here it scales with the log-time ~ A _1 | logfi.| 3> 
A -1 . This makes the coherent test state a factor of 
T X/y/n ~ | log % | times better than any of the single 
Gaussian test states. 

We can also look at the spectral envelope S^ n which 
is the fourier transform of the short-time autocorrelation 
function (see Eq. |l6| ): 



JEt 



t 



(59) 
(60) 



where again in the last line the limiting (% — > 0) form 
Eq. [55] has been used. The energy envelope is centered 
at E = (j>, just like the smoothed single-wavepacket local 
density of states, but the the peak is both narrower and 
taller by a factor scaling as | log7i|. 



B. Numerical Tests 

We now check the results obtained in this section, using 
again the ensemble of kicked-baker systems introduced in 
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Section III C . The short periodic orbit will again be the 
fixed point of the middle baker's strip, with exponent A 
set by the logarithm of the width of this strip. We begin 
by looking at the smoothed local densities of states 5* n 
and S\ in , for the universal test state and the simple 
Gaussian \a), respectively. The width a of the starting 
Gaussian is set to y/w-y, so that the aspect ratio of the 
Gaussian is equal to the aspect ratio of the rectangular 
middle strip in which the classical dynamics is lineariz- 
able. The wavepacket can then expand the same number 
of steps in either time direction before reaching the edge 
of the linear regime. The test state *f? is constructed using 
a cutoff set by c = 0.6 (see Eq. 

Ensemble averaging is performed over many kicked- 
baker systems of the same (reasonably large) exponent 
A = | log 0.18| , and of system size (Hilbert space dimen- 
sion) TV — 1/h — 800. Local densities of states for \a a ) 
and | "J) are ensemble-averaged and smoothed, with the 
resulting envelopes plotted in Fig. ||. Theoretical curves 
obtained form Eqs. [ll], ^9] are also plotted for compari- 
son. Excellent agreement is observed between the data 
and the predictions based on the linearized theory. Fur- 
thermore, we see that the spectral envelope for the hyper- 
bolic test state 'J is significantly narrower and taller than 
the corresponding envelope for the Gaussian wavepacket, 
again in accordance with prediction. We should note here 
that the hyperbolic test state is constructed here with the 
very modest log-time cutoff To = 0.90. There are three 
reasons for the smallness of To in this example: 1) the 
stretching factor e A « 5.6 is rather large, 2) the system 
size, and particularly the size of the linearizable region, 
are modest, 3) and finally the free parameter c has been 
set at a rather conservative (small) value. With regard 
to the last point, we should note in particular that in- 
creasing the cutoff parameter c (Eq. |24|) will give rise to 
a sharper envelope, with larger inverse participation ra- 
tio, though at some cost to the accuracy of the formulas 
Eqs. |57], |5^, etc. In effect, there is a tradeoff between 
keeping the test state well inside the linear region and 
thus being able to obtain with good accuracy its statis- 
tical properties (smaller c) versus maximizing the local- 
ization properties of the hyperbolic test state by allowing 
it to some extent to leak out of the linear region (larger 
c) . All of this will become clearer as we go on to discuss 
IPR measures for the universal test states. Of course, 
none of these 0(1) considerations affect the basic scaling 
predictions: namely the height, inverse width, and IPR 
of the spectral envelope for <3> all scale inversely with A 
for small A and also logarithmically with 1/h for small 
h. In particular, the hyperbolic (coherent 4" test state) 
spectral envelope gets arbitrarily taller and narrower as 
h — > for a fixed classical system, while the correspond- 
ing spectral envelope for a single Gaussian packet test 
state remains unchanged. 

We next probe the behavior of the mean IPR^ as 
a function of TV and c (with exponent A again set to 
| log0.18|). In Fig. |, we plot the IPR versus cutoff time 
T , for five sets of data: TV = 1/h = 50, 100, 200, 400, and 



800 from bottom to top. For each value of TV, 26 values 
are plotted: from left to right c = (1.1) J , j = —20 . . . + 5. 
The five values at T = represent simple Gaussian test 
states (c — ► 0). The upper dashed curve represents the 
theoretical prediction of Eq. |57], which should hold for 
large values of TV. Good agreement with the data is 
obtained for TV > 200, as finite-size effects become less 
relevant The six rightmost points on each data curve 
represent c > 1, and some deviation from the linear the- 
ory prediction is expected to start setting in there. We 
also note that at very large values of To (requiring cor- 
respondingly larger values of log TV), the theoretical pre- 
diction converges to the linear asymptotic form of Eq. |58| 
(lower dashed line). 

We note for purposes of comparison that the single- 
wavepacket scarring strength IPR a is predicted to be 
3.66 for this value of A (see Eq. |2(], noting that a quan- 
tum fluctuation factor F = 2 is appropriate for complex 
eigenstates). This indeed is close to the value attained by 
the single wavepackets (T = 0), at least for TV > 200. [] 
We see that IPR values significantly larger than this can 
be attained using the hyperbolic test states, especially 
for larger values of TV. The data is consistent with the 
prediction that the IPR (at fixed value of c ~ 1) scales 
logarithmically with TV for large system size TV. 

Finally, to close this section we present in Fig. || the 
Husimi phase-space plot for the hyperbolic test state \F 
In this figure, To has been taken to be very large com- 
pared to A^ 1 , i.e. the linearizable regime is much larger 
than a unit Planck cell, and also much larger than the 
phase-space area shown in the figure. We note that the 
phase-space Husimi picture is universal, and in particu- 
lar independent of the exponent A, since a change in the 
value of A in Eq. |5l| can of course always be absorbed 
into a redefinition of To and the overall normalization TV. 
In other words, in the log 4 — * oo limit, the state \f de- 
pends only on the linear region size parameter A/h, and 
when we further look well inside the area A, its structure 
is completely free of any parameters. The phase space 
area shown in the figure is i.e. it contains 

144/27T Planck-sized cells. We see from the Figure that 
^ lives at the periodic point in the center of the plot (the 
size of the bright region at the periodic point being set 
by h) , and symmetrically along the linearized stable and 
unstable manifolds. This picture will be important to 
us for comparison purposes when we study off-resonance 
universal test states in the following section. The ana- 
lytic expression used to obtain the density plot in Fig. ^ 



4 Notice that because the width of the central strip here is 
quite small (0.18), TV = 200 corresponds to a size of only 
0.18 x 200 = 36 for the linearizable region (in units of h). 

The deficit in the measured single Gaussian IPR values 
(To = 0) for finite TV, i.e. the extent to which these fall 
below the limiting value of 3.66, is in close correspondence 
with similar deficits in the hyperbolic state IPR's (To > 0). 
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will also be given there (Eq. 



V. OFF-RESONANCE AND OFF-ORBIT 
SCARRING 

Going back to the construction of the hyperbolic test 
state \& in Eq. pi, we notice that there the Gaussian 
wavepackets are all added in phase, giving rise to a pre- 
ferred energy <f> which is the same as that for any single 
wavepacket considered individually. We may, however, 
equally well consider the more general form 



* = A/" dte 



"* ' T Se m \a n 



(61) 



where 9 is an arbitrary phase accumulated per time step. 
This extra phase should give rise to a state that prefers 
to live at an energy different from the one that exactly 
quantizes the periodic orbit (i.e. E = <ft). In turn, this 
energy shift may be expected to give rise to phase space 
structures that lie away from the invariant manifolds of 
the periodic orbit, i.e. above and below the separatrix 
constructed from these manifolds j2j|. These intuitive 
expectations turn out to be justified, as we now shall see. 

We begin once again by computing the short-time (lin- 
earized) autocorrelation function. The generalization of 
Eqs. pil, |55| works out to be: 



A* n (t) = Qe' l4>t 



dy 



e -(t-if/Tg-ie(t-%) 



Vcoshy 



-i{4>+e)t p -t 2 /T* 



(62) 
(63) 



In the last line, the limits A — * and To A ~ lo g 4 — > 00 
have been taken, in complete analogy with Eq. p5| . 

^From these expressions, we proceed to obtain spectral 
envelopes and IPR values, exactly as we did previously 
for the 9 = special case. It is interesting to note here 
that in the asymptotic limit of Eq. the spectral en- 
velope will be exactly the same as that obtained previ- 
ously (Eq. ^) , but shifted so as to be centered at energy 
E = <fi + 9. Because the shape of the envelope is un- 
affected by the value of 9 in this limit, the IPR is still 
asymptotically given by our previous formula, Eq. It 
is important to note, however, that for realistic values of 
log 4 * this asymptotic regime may not be attained, and 
the more general formulas Eqs. |57|, [39], |6^ should be used 
instead. 

In Fig. ^, we plot (numerically obtained) smoothed 
spectral envelopes for the hyperbolic states V^, for the 
same ensemble of systems as was used in producing 
Fig. ||. The envelopes correspond to 9 = (the tallest 
envelope, already seen previously in Fig. ||), through 
9 = —2ir, moving to the left in steps of 27r/20. Even 
though the asymptotic form of Eq. ^33] predicts all the en- 
velopes should have the same shape, being merely shifted 
to the left by angle 9, in reality we see this is not quite 



the case for finite values of \og(A/h). Because of the fi- 
nite linearizable volume of phase space, the off-resonance 
hyperbolic states have significantly less well-defined en- 
velopes compared to the 9 = state. That is because 
the asymptotic form assumes most of the autocorrela- 
tion function comes from long-time overlaps of wavepack- 
ets with very different widths. At finite system sizes, a 
very important correction is the partial self-cancellation 
in W coming from wavepackets of comparable widths 
being added together with very different phases. This 
correction is, of course, taken into account in the more 
general form of Eq. 



62 



which does in fact predict less 
sharp envelopes (and consequently lower IPR's) for the 
off-resonance states. The important point to notice here, 
however, is the presence of a very significant localiza- 
tion effect even for \9\ > A, i.e. at energies well out- 
side the resonance of the original Gaussian wavepacket. 
States at such energies are not strongly scarred accord- 
ing to the original (Gaussian wavepacket) definition, nor 
are they particularly enhanced along the stable and un- 
stable manifolds of the orbit, as measured by the on- 
resonance hyperbolic test state. However, such states 
do have enhanced intensity relative to the off-resonance 
hyperbolic states ^(9 ^ 0), which live in hyperbolic re- 
gions on cither side of the separatrix (see Fig. |l^ later 
in this section). Although the enhancement factors for 
such states are quite modest for the parameters chosen 
(~ 2 over RMT), they will of course grow with A -1 and 
with log A/h, as we discussed earlier. In particular, in 
the high energy limit (h — > 00) of a given classical sys- 
tem, scarring by the off-resonance states is expected to 
be equally strong compared to the on-resonance (9 = 0) 
form of scarring. 

To demonstrate the preceding assertion, we plot in 
Fig. |^ the (theoretically computed) IPR as function of T , 
for several values of the off-resonance angle 8. To facili- 
tate comparison with Fig. ||, we again choose stretching 
exponent A = | log 0.18| . The six curves represent (from 
top to bottom at To = 1) values of 9 from to 7T, in steps 
of 7r/5. The top curve has thus already appeared previ- 
ously in Fig. ||. The lowest curve will later be compared 
with data in Fig. ||. The asymptotic form of Eq. [38], to 
which all these curves converge at large T (large | log K\), 
is shown in the figure as a dashed line. 

An interesting point to notice here is that for suffi- 
ciently large 9, the IPR can drop at moderate Tq from 
its single- wavepacket value (at To = 0), before eventu- 
ally recovering at larger To (this is a result of the self- 
cancellation effect alluded to earlier). However, for (ex- 
ponentially) large systems, scarring is equally strong for 
the different values of the off-resonance parameter 9. In- 
tuitively, for a large system (or small h), the size of 
the linearizable region in which the test state W is con- 
structed is very large compared to the size (~ h) of a 
single wavepacket. Thus, a given wavepacket used into 
the construction of \& in Eq. [H] has very little overlap 
with most of the other wavepackets in the sum. For 
this reason, the cancellation effect arising from similar 
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wavepackets being added together with different phases 
(for non-zero 9) becomes insignificant in the semiclassical 
limit. 

To make contact once again with the data, we choose 
9 = tt, and in Fig. @ do for this case what we did for 
9 = earlier in Fig. |f Again, data for N = 50, 100, 200, 
400, and 800 is plotted, and the same range of c-values is 
used. Only for N > 200 do we see the recovery towards 
larger IPR values that we expect ^from the theoretical 
curve (dashed). Recall that 3.66 is the single- wavepacket 
IPR value. 

In Fig. ||, the system size N = 200 is fixed, with cutoff 
parameter c = 0.6, and the IPR is plotted versus the off- 
resonance angle 9, ranging here from to 2tt. Agreement 
between theory and data is in this case (surprisingly) 
good. We see for this value of TV that by 9 — 2ir, the 
IPR values are already approaching the RMT value of 2, 
indicating almost no localization. 

Finally, we return to the Husimi representation of \f r , 
which we began to discuss already in the previous sec- 
tion. Husimi representations of the off center scarring 
were discussed in p{| , where the inverted oscillator eigen- 
states were probed with Husimi states. Here, we have 
shown how to generate test states which are sensitive to 
off-center scars developing in chaotic systems. For ref- 
erence, the expression describing the Husimi intensities 
is: 
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(64) 



Here 



go, pa 



is a phase-space Gaussian centered at phase- 



space point (qo,Po) of width yh in both position and 
momentum: 

g q0tP0 (q) = (7r7i)-V4 e -(g-go) 2 /2ft+ipo(9-<2o)/ft _ (6S) 

Recall that we are working in a coordinate system where 
the fixed point is located at the origin (q,p) = (0,0), 
and the linearized invariant manifolds are horizontal and 
vertical (Eq. ||) . As the linearizable region becomes large 
compared to the phase space region of interest (To — > 
oo), the Husimi density of Eq. ^ depends on only three 
dimcnsionlcss numbers: an off-resonance parameter 9/X 
and the phase space coordinates qo/Vn~, po/Vh. 

Previously, in Fig. ^, we have seen the Husimi density 
of ^ for 9/X = plotted in the square phase space area 
—6 < qo/y/ft^po/Vfo < +6. In Fig. [l(] we present the 
analogous picture for (a) 9/X = 0.8, and (b) 9/X = 2.5. 
As the energy goes further off resonance, the state \E' 
moves away from the periodic orbit and its invariant 
manifolds, and shifts into two of the quadrants separated 
by the manifolds (qp > or qp < 0, depending on the 
sign of 9). For large \9\/X, narrow hyperbolic regions 
in phase space are accessed, lying further and further 
from the periodic orbit itself. For 9/X ~ ±1, 'J lives 



near the hyperbolas qp ~ ±h surrounding the periodic 
point. This regime corresponds to a spectral envelope 
for ^ which is centered at energy A away from the EBK 
energy, i.e. at the edge of the single-Gaussian spectral 
envelope. Eigenstates having strong overlaps with 4" are 
now barely scarred at the periodic point itself. If we go 
further into the regime \9\/X ^S> 1, the spectral envelope 
of ^ is now centered at an energy which is outside the 
envelope of the single wavepacket. Then, states overlap- 
ping such a test state $ will have stronger than expected 
intensity on hyperbolic regions surrounding the periodic 
orbit, but will not be scarred at all (and may even be 
antiscarred) on the orbit itself. 

The connection discussed above between an energy 
shift away from the EBK value and hyperbolic phase- 
space structures may be understood very simply by con- 
sidering the evolution of off-center Gaussians g qo , Pa of 
Eq. |6^. Such a Gaussian is not an optimal test state for 
measuring scarring, because the autocorrelation function 
decays quite rapidly, especially for q$ + p\ ^> h (to be 
contrasted with 'J, which lives along the entire hyper- 
bolic region, and thus has much larger self-overlaps at 
short times). However, the phase information in the au- 
tocorrelation function for <? go ,p is quite relevant: 



A hn (t) = e 



, ,f "n+ p o ""h 2 At/2 
— p LC P L p h cosh At 



(compare with Eq. ^). Upon fourier transforming this 
to obtain a spectrum, we obtain an expression for the 
optimal energy as a function of phase-space location: 



E- <t> w A 



qp 



(67) 



Thus the off-resonance parameter 9 of our test state \T/ is 
then expected to be associated with phase-space hyper- 
bolas: 



qp 

h 



(68) 



VI. EXTENSION TO HIGHER-PERIOD ORBITS 
AND CONTINUOUS TIME 

A. Longer orbits of maps 

The analysis of the previous two sections has focused 
on fixed-point periodic orbits, but it generalizes in a 
straightforward way to longer orbits and to continuous 
time. We shall see below that the benefits of using uni- 
versal scar measures instead of simple Gaussians become 
even greater when longer orbits are considered. 

Consider again a periodic orbit (of a map) of period 
Tp, with periodic points x p (p = 0...Tp — 1), as in 
Eq. ^6[ Let —0 be the phase accumulated over one full 
iteration of the orbit, and A the corresponding stability 
exponent. The short-time autocorrelation function for a 
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Gaussian \a) centered at any of the periodic points is a 
generalization of Eq. ^: 



A na (t) 
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(69) 



The IPR as function of A is then the same as for a 
period-one orbit (because it is given simply by the sum 
of the short-time return probabilities, Eq. while the 
short-time spectral envelope Su n has Tp peaks in the 
quasienergy domain [0, 2ir], of height scaling as A^ 1 , and 
width scaling as X/Tp. The peak energies are of course 
those that 'quantize' the orbit: TpE = 0mod27r, or 



E k = 



2nk 



, k = . . . T P - 1 



(70) 



We notice that both the maximum scarring strength and 
the IPR can be large only for small A, which becomes 
difficult to achieve for the longer orbits Tp > 1 (normally 
A grows linearly with Tp). 

We now proceed to construct the universal test state 
iff for such an orbit, having made a choice of quantization 
energy E k : 



xei( WpWP^| aifrteAt/Tj>) . (71) 

Here f p is a stretching factor, and 4> p is a phase, both 
defined by 



c —i<t>p 



(72) 



f p and <fip take into account the fact that stretching and 
phase accumulation along the orbit may both be non- 
uniform; of course, fr P — 1 and 4>t p — 4>- The factors 
fp are of order unity and thus not very important in the 
semiclassical limit log A/h — ► oo when the linearizable 
region is very large; on the other hand, the phases <p p are 
crucial for getting constructive interference. The param- 
eter To is defined as before (Eq. ^J) , using the area of the 
linearizable region around periodic point xq, and 9/X is 
an off-resonance parameter, as discussed previously. 

The short-time autocorrelation function of the state iff 
then has the same form as what we found previously for 
the special case Tp = 1 (Eqs. H E35], E32], E33|), replacing 



To -» T Tp , A - X/Tp , 9 -> 9/T P 



E k (73) 



throughout. Note that A and 9 are defined as stretching 
exponent and phase per orbit period rather than per time 
step, and likewise To is the log-time measured in units of 
the orbit period. Eqs. [59], |6(] describing the shape of 
the linear spectral envelope, and Eqs. |57], |58| for the IPR 
of the universal test state iff undergo the same simple 
modifications and are then applicable to the case of Tp > 
1. 



Let us compare these results with ordinary (single 
Gaussian) scarring for Tp > 1 as well with the scarring 
of a fixed point, which we focused on in the previous two 
sections. 

For general Tp > 1, the autocorrelation function of a 
Gaussian wavepacket \a a ) has 0(A _1 ) strong recurrences 
spaced Tp steps apart, and thus stretching over a total 
time scale of order TpA -1 . In the quasienergy domain, 
this leads to a set of Tp evenly spaced spectral envelopes, 
each with width, height, and IPR scaling (for small A) as: 
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(74) 



Although the width scales with the exponent per time 
step, due to the presence of Tp of these envelopes, the 
maximum expected scarring strength and the IPR both 
scale only with the total exponent per iteration of the 
entire orbit, and thus are expected to deviate less and less 
from perfect ergodicity as longer orbits are considered. 

Let us repeat the same analysis for the state iff, which 
takes properly into account all of the linearized dynam- 
ics around the periodic orbit. The short-time autocor- 
relation function does not decay until a time of order 
T)Tp ~ TpA -1 log A/h. This produces a single peak cen- 
tered at quasienergy E k (or shifted by an off-resonance 
phase 9/Tp). The width, height, and IPR scale as 



IPR- 



w a ~ X/{T P log A/h) 
' h a ~ (Tp log A/h)/ X 



(1 envelope, centered at some E k + — ) . (75) 

Tp 

Apart from the logarithmic enhancement which leads to 
more and more deviation from RMT in the h — > limit, 
we also notice that only one peak is present (a choice 
of quantization energy E k having been made), and all 
quantities now depend only on the ratio X/Tp, the sta- 
bility exponent per time step. This measure of scarring 
therefore allows us to see strong effects even for longer 
periodic orbits, as long as the stretching per time step 
remains moderate. 

If the total exponent A (and not just X/Tp) is small, we 
can equivalently use the linearized dynamics to generate 
our universal test state (compare Eqs. ||(], |5^). Eq. |52| 
generalizes easily to 



dyn 



t 



z -t 2 /T P T$ ei (E k+ e/T P )t laiin{t)h (?6) 



where \a) is a Gaussian wavepacket of width a centered at 
any point along the periodic orbit. Eq. |7l] can be thought 
of as an averaging of Eq. fr| over the initial width a from 
some Co to ooe A . For small A, the averaging procedure 
is unnecessary, all the states being essentially identical, 
and the simple expression of Eq. [76] well describes the 
universal test state for any choice of a. 
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B. Hamiltonian systems 

The entire analysis can be applied also to Hamiltonian 
systems in an essentially unchanged form. Let Tp again 
be the period of the orbit (now measured in real time 
units rather than in time steps), and let A and <fi still be 
the exponent and phase per one iteration of the orbit. 
Then an optimal Gaussian centered anywhere on the or- 
bit and aligned along the stable and unstable manifolds 
(with width a along the unstable manifold) has a linear 
autocorrelation function given by a sum over iterations 
of the orbit: 



A nn (t)=J2w(t-nT P ) 



V cosh Ar 



(77) 



The very-short-time window function w(t) describes the 
self-overlap of the wavepacket as it intersects itself once 
every period; it is associated with the nonzero width a x 
of the wavepacket in the direction of the orbit, and has 
a scale At ~ <J x /v ~ a x /\[E. (In the last equality, the 
particle mass has been assumed to be unity, as it will 
be throughout.) Fourier transforming, we obtain a very 
wide envelope (of width AE ~ %y/E/a x ) associated with 
the energy uncertainty of the wavepacket itself. Multi- 
plying this are the scarring envelopes, centered at the 
quantizing energies, of widtl^j 



and separated by 



W a ~ hX/Tp 



h/Tp 



(78) 



(79) 



The normalized peak height, and the IPR are therefore 
given by 



IPR ~ h a ~ A" 



(80) 



Again we see that ordinary measures of scarring are typi- 
cally unable to resolve scarring arising from longer orbits, 
because the exponent A for such orbits is generally not 
small. 

The universal state ^ is constructed by analogy with 
Eq. fn] as 



■^=jvj^ dr J dte- t2/T P T » 



i(E k +he/T P )r/h+iet/T P -i4, Tl ^ 



51) 



Here r is a time parameter parametrizing the orbit x T , 
f T and 4> r are as before a stretching factor and phase 



Notice that E is now a real energy, rather than the Floquet 
phase it was for a map, and Tp has units of time instead of 
step number, hence the factor of h in the equations following. 



associated with short time evolution from xq to x T , Eu 
is a quantization energy, and 9/X is again an optional 
off-resonance energy shift parameter. For small A, this 
can be written as 



dyn 



E 



-e/Tlrl e i { E k+ ne/T P )t/n^ t)) (82) 



The short-time autocorrelation function again has a 
decay time scale TpTo ~ TpA -1 log A/% (with no window 
w(t) present), leading to a single spectral envelope cen- 
tered at energy E/. (Ek + %9/Tp for off-resonance states), 
and of width and height given by 



wq, ~ %\/{T P log A/%) 
hy ~ {T P log A/%)/%\. 



(83) 



The IPR is somewhat difficult to talk about in this case 
because IPR's (like any other measure of quantum er- 
godicity) can only be measured relative to some already 
known energy window which takes into account various 
conserved quantities. In this case, the only plausible 
window is the spectral envelope of the original Gaussian 
wavepacket of width a x in the direction of the periodic 
orbit (see discussion following Eq. |7?]). Then the IPR is 
given by 



IPR* 



ET P j A 
Xa r ^ % 



(84) 



Notice that the combination X/yETp is just the expo- 
nent per unit length of the orbit. Of course, our result 
depends on the width a x of each Gaussian along the or- 
bit. The enhancement of the IPR for \& over the corre- 
sponding IPR for the single Gaussian (Eq. |8(]) can also 
be thought of as being given by the usual logarithmic fac- 
tor times the ratio y/ETp / a x of the length of the orbit 
compared to the wavepacket size. 

In any case, we again see that long orbits can be easily 
resolved using this improved scarring measure. A suffi- 
cient condition to get significant enhancement is for the 
exponential stretching to be small on the time scale of 
the very-short-time window w(t), i.e. during the time 
it takes for the wavepacket to traverse its width while 
moving along the orbit. This criterion is of course inde- 
pendent of the orbit length. Even if the criterion above is 
not satisfied at low energies, one nevertheless gets strong 
scarring in the semiclassical limit. For a two-dimensional 
billiard system, the increase in scarring strength (as mea- 
sured by peak height of the spectrum or by the IPR) 
scales with energy as 



\f~E\oeE 
IPR* - E lfi \ogE. 



(85) 
(86) 



(The power-law IPR scaling arises ^from the wavepacket 
width a x scaling with energy as u x ~ E^ 1 / 4 . This is a 
natural scaling which keeps the a;— uncertainty and the 
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uncertainty in the x— momentum p x in a fixed ratio rela- 
tive to the total size of the accessible phase space. With- 
out this scaling the increase in IPR with energy would 
be only logarithmic.) 



long the orbit period remains small compared to the log- 
time. Orbits of arbitrary length can therefore be strongly 
scarred using this measure, provided a correspondingly 
small value of fi is chosen. 



C. Numerical tests 



VII. SUMMARY AND CONCLUDING REMARKS 



We conclude this section with a numerical example of 
the localization enhancement obtainable for longer orbits 
using the universal test-state approach. For this purpose, 
we choose a modified versio n of the kicked-baker system 
introduced in Section QIC. Instead of having the kicks 
act on the left and right strips of the three-strip system, 
we have one act now only on the middle strip, leaving 
the left and right strips to undergo ordinary bakcr-likc 
dynamics, i.e. horizontal expansion and vertical com- 
pression. Any periodic orbits contained entirely in these 
two side strips thus have locally orthogonal stable and 
unstable manifolds (of the form Eq. ||), with a stretch- 
ing exponent and action given simply in terms of the 
widths wq and W2 of the left and right strips. In par- 
ticular, consider the period-2 orbit that jumps from the 
left strip to the right strip and back. Its periodic point 
in the left strip is given by q = wq(1 — W2)/(l — wo^), 
p = (1 — W2)/{1 — wqW2)- The other periodic point is 
obtained by interchanging the q and p coordinates. The 
stretching exponent for one full iteration of this orbit is 
A = | log wq wi |, and the corresponding phase is given by 
4> = u?o(l — W2) /(l — WoW2)h. Thus, a desired value 
for the exponent and phase can be fixed by selecting the 
three baker strip widths, and the kick strength acting on 
the middle strip is then used to provide ensemble aver- 
aging over the details of the long orbits (nonlinear recur- 
rences) . 

We select for our example widths u> = 0.40, u>2 = 0.42, 
as before, leading to an exponent A = | log 0. 168| for our 
chosen orbit, and work with the matrix size N = 800. In 
Fig. [ll], the smoothed local density of states for a Gaus- 
sian wavepacket, obtained by averaging over several real- 
izations, is represented by the double-peaked solid curve. 
The theoretical prediction, given by the linearized dy- 
namics of Eq. ^39], is shown by the dashed curve. The nar- 
row, single-peaked solid curve centered at one of the two 
possible quantization energies is the similarly smoothed 
local density of states for the universal state con- 
structed once again with c = 0.6 and 6 — 0, as in Fig. [|. 
Again, the corresponding dashed curve is the theoretical 
prediction based on the linearized dynamics and agrees 
well with the data. Notice that in this case the difference 
between Gaussian and universal scarmometers is more 
dramatic than in Fig. ^, the reason being that scarring as 
measured by '5 depends only on the stretching rate along 
the unstable manifold per unit time, not per iteration of 
the entire orbit. While scarring strength as measured 
by a single wavepacket drops off with orbit length, scar- 
ring strength as measured by is length-independent as 



In this paper we began by reviewing the theory and ex- 
isting measures of scarring. We were then able to estab- 
lish considerably refined and arguably universal scarring 
measures. The refinement means in practice that much 
larger deviations from RMT behavior are predicted us- 
ing the refined test states. The test states are not special 
from the standpoint of random wavefunctions, but they 
pick up structures which exist in eigenstates of dynamical 
systems. 

The two major issues of generalization of scar mea- 
sures which we faced are 1) coherence (or lack of it) 
in superpositions of localized wavepacket states, and 2) 
summing over all the points of periodic orbits whose pe- 
riod is greater than one iteration of a map. The smooth 
Hamiltonian version of this is to coherently add up pack- 
ets all along the orbit, making phase space tubes which 
are related to the tubes of Agam and Fishman p2fl ; see 
also Uj. 

The universality mentioned above stems ultimately 
from the use of the linearizable domain near periodic 
points in the construction of scar measures. The scar 
test states are the optimal ones which can be constructed 
with the linearized dynamics. In turn, we argued that 
the linearizable portion of the dynamics was a reason- 
able stopping point for the definition of scar strength. 
Going beyond the linearizable dynamics is certainly pos- 
sible and semiclassically viable, but a problem arises in 
that one begins to approach the construction of individ- 
ual eigenstates, at least in favorable cases ]l8| |, which is a 
somewhat disturbing limit. The reason this is disturbing 
is that such "test" states for scarring pick up (in the ideal 
limit) only one state, which brands the whole eigenstate 
as a scar. Moreover, pieces of classical manifolds far from 
any given periodic orbit will have been incorporated in 
the longer time dynamics of such a test state. Indeed 
it is not at all clear that any one periodic orbit should 
dominate the others in such a state. These new periodic 
orbits would begin to play a role in the long time dynam- 
ics (on the order of the Heisenberg time), so we would 
not even be speaking of a scar of a given periodic orbit. 
Given all these factors it seems reasonable to stop at the 
linearizable zone surrounding given periodic orbits. 

Various numerical tests made possible by ensemble av- 
eraged baker map results supported the measures estab- 
lished here. The enhancements in IPR made possible by 
our optimized coherent measures can be modest (factors 
of 1.5 or 2) for reasonable fi and short orbits, but much 
more dramatic for longer period orbits, as compared to 
a single Gaussian wavepacket measure. 
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Finally, we have given a theoretical basis for the "off 
center" scars living on the hyperbolic manifolds near (but 
not on) a given periodic orbit, and provided them with 
test states sensitive to their presence. 
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FIG. 1. The correlation C ab (Eq. ||) between local densi- 
ties of states for two wavepackets \a) and \ b), lying on the same 
periodic orbit of instability exponent A, is plotted as a func- 
tion of width parameter z. The upper and lower curves cor- 
respond to A = log 5 and A = log 10, respectively. The width 
(along the unstable manifold) of wavepacket \b) is e A ' n+z ' 
times that of \a), where n is an integer. Numerical data 
for t he en semble of kicked-baker systems described in Sec- 
is plotted along with the theoretical curves from 



tion 



HIC 



Eq. g^! As the stretching factor e A gets closer to unity, so 



does the correlation C a b 



FIG. 2. The minimum correlation in the local densities 
of states Cab for two wavepackets on an orbit of exponent A 
is plotted versus the scarring IPR enhancement factor g{\) 
for such an orbit. The theoretical curve is obtained from 
Eqs. wh ile the data again comes from the ensemble 

of Section ElIC. The minimum inter- wavepacket correlation 



gets very close to unity for significant scarring enhancement 
factors g(X). 



FIG. 3. Smoothed local densities of states are plotted 
for the universal hyperbolic test state !<&} (higher peak) and 
a Gaussian wavepacket \a) (lower peak), on a periodic orbit 
with exponent A = | log 0. 18| . The system size is N = 800. 
Cut-off constant c = 0.6 (see Eq. ^ij) is used to construct the 
state |^}. The theoretical curves (dashed) are obtained by 
fourier transforming the linearized autocorrelation functions 
of Eqs. [^J |^, respectively, while the data (solid curves) is 
obtained by ensemble averaging. 



19 



FIG. 4. In this figure, the inverse participation ratio (IPR) 
for hyperbolic test state |$) is plotted versus the log-time 
cutoff To (see Eqs. [HI], for various values of system size 
TV. From bottom to top, the five curves correspond to TV = 50, 
100, 200, 400, and 800. For each TV, 26 points are plotted, 
for c = (1.1)"', j ~ —20... + 5. The orbit has exponent 
A = | log0.18|, as in the previous figure. The upper dashed 
curve is the TV — > oo theoretical prediction (Eq. |57], F = 2), 
which converges to the asymptotic prediction of Eq. ^ (lower 
dashed line) for large To. The linearized theory is expected to 
start breaking down for c > 1 (rightmost six points on each 
data curve). 



FIG. 5. This figure is a Husimi plot in phase space of 
the universal hyperbolic state | ^} for log-time cutoff T — > oo 
(i.e. the region plotted is well inside the linearizable region in 
classical phase space). The fixed point it at the center of the 
plot, and the horizontal and vertical axes correspond to the 
unstable and stable manifolds, respectively. The total area of 
the plot is 12y/K x 12y/h. 



FIG. 6. Smoothed local densities of states are plotted 
for the off-resonance universal hyperbolic test state for 
off-resonance angle 9 ranging from (tallest peak) in steps 
of 7r/10 (to the left), through — 2%, on a periodic orbit with 
exponent A = | log . 1 8 1 . As in Fig. ^, the system size is 
TV = 800, and the cut-off constant c is set to 0.6 (see Eq. [mJ). 
The data is obtained by ensemble averaging. 



FIG. 7. Theoretically computed IPR values for the uni- 
versal hyperbolic test states ^ are plotted versus the log-time 
To for several values of off-resonance angle 9. From top to 
bottom at To = 1, the solid curves represent 9 — 0. . .71", in 
steps of 7r/5. The dashed line is the limiting value for all 
of these at large To (Eq. |H^ ). The single-Gaussian IPR (the 
To — * limit) is 3.66 for this value of A. 



FIG. 9. Here the IPR is plotted versus off-resonance an- 
gle 9, for A = | log 0. 18| , TV = 200 and cutoff time param- 
eter c = 0.6. The solid curve with errorbars represents en- 
semble-averaged data, while the dashed curve is the theory 
(Eqs. |62| ). At small angles, the scarring strength is signif- 
icantly above the single wavepacket IPR value of 3.5 (theo- 
retical prediction: 3.66), while for larger angles we approach 
the ergodic RMT value of 2. 



FIG. 10. Analogous to Fig. jsj this figure presents Husimi 
plots of the universal test state ^ for off-resonance parameter 
values (a) 9/X = 0.8 and (b) 9/X = 2.5. As before, the lin- 
earizable region is taken to be much larger than the displayed 
area of size 12\/ft x 12y/h, centered on the periodic orbit. 



FIG. 11. Smoothed local densities of states are shown 
here for a Gaussian wavepacket placed on a period-2 or- 
bit (double-peaked solid curve), and the universal test state 
constructed on the same orbit (tall single peak). The 
dashed curves represent theoretical predictions based on the 
linearized dynamics near the periodic orbit in question. The 
system is a kicked baker's ma p with kick potential acting on 
the middle strip (see Section VI C ) , and the periodic orbit 
has a total exponent A = | log 0.168| over the two-step period. 
One of two possible on-resonance energies has been chosen 
for the test state which is again constructed using cutoff 
constant c = 0.6 (as in the analogous calculation in Fig. ^ for 
a period-one orbit). The enhancement here is more dramatic 
due to the fact that universal scarring strength depends only 
on the exponent per unit time along the orbit, not on the 
orbit length itself. 



FIG. 8. The IPR is plotted here as a function of log-time 
To, as in Fig. ^, but for the off-resonance test state "^(9 — ty) 
of Eq. [n]. Again, five data curves corresponding (bottom to 
top) to TV = 50, TV = 100, TV = 200, TV = 400, and TV = 800 
are shown in the Figure, with cutoff parameter c varying from 
(l.l)" 20 to (1.1) +5 from left to right within each curve. The 
stability exponent is A = | log 0. 18| , as before. The upper 
dashed curve is the TV — > oo theoretical prediction obtained 
from Eqs. |57|, p2| . The lower dashed line is the asymptotic 
form of Eq7]58]7to which the upper curve converges in the 
To — > oo limit. The linearized theory is again expected to 
break down for c > 1 (the rightmost six points on each data 
curve) . 
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